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Non-thermal Electrons at the Earth's Bow Shock: A 'Gradual' Event 
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Earth's bow shock is known to produce non-thermal electrons which are generally observed as a 'spike' in 
their flux profile. Here, in this paper, we present an analysis of electron and whistler wave properties for a quasi- 
perpendicular shock crossing that is supercritical, but subcritical to the so-called whistler critical Mach number. 



above which whistler waves cannot propagate upstream. We have found that the amplitudes of whistler 



waves increased exponentially as a function of time prior to the shock encounter, while the suprathermal (> 2 
keV) electron flux similarly increased with time, although with differing e-folding time scales. Comparison of 
the electron energy spectrum measured within the ramp with predictions from diffusive shock acceleration theory 
was poor, but the variation of pitch angle distribution showed scattering of non-thermal electrons in the upstream 
region. While not finding a specific mechanism to account for the electron diffusion, we suggest that the whistlers 
seen probably account for the differences observed between this 'gradual' event and the 'spike' events seen at 
shocks with no upstream whistlers. 
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1. Introduction 

Energetic electrons with energies more than 20 keV have 
been observed at and near the Earth's bow shock (e.g. Fan et 
al., 1964; Frank and Van AUen, 1964; Anderson, 1965, 1969; 
Vandas, 1989). Since larger electron flux can be found on the 
interplanetary magnetic field (IMF) tangent to the bow shock 
(e.g. Anderson et al., 1979; Kasaba et al., 2000), electrons 
are considered to be accelerated in the quasi-perpendicular 
region where the shock angle Obu is larger than 45°. Gosling 
et al. (1989) was the first to carry out comprehensive anal- 
ysis of suprathermal (<20 keV) electrons across the shock 
front. In their quasi-perpendicular shock events, energetic 
electron flux was enhanced at the shock transition, and be- 
cause of the localized feature, they termed their events as 
'spike' events. The energy spectrum showed a power-law 
form with the spectrum index of 3-4. The pitch angle dis- 
tribution was almost isotropic at the transition layer while 
it was anisotropic in both the upstream and the downstream 
regions. More recently. Oka et al. (2006) conducted a statis- 
tical analysis of the power law indices measured in the shock 
layers. They reported that the power-law index of electron 
energy spectra is regulated by the so-called whistler critical 
Mach number M™;^, which is defined as the critical point 
above which whistler waves cannot propagate upstream. 

In this paper, we report a shock crossing event that showed 
'gradual' profile of non-thermal electron flux in association 
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with intensification of precursor whistlers. Contrary to the 
spike events reported by Gosling et al. (1989), the electron 
flux increased exponentially with decreasing distance from 
the shock. The event has been determined to be subcritical 
in relation to the whistler critical Mach number M™;^. We 
will describe properties of the waves and discuss origin and 
transport of the non-thermal electrons in this gradual event. 



2. Observation 

Our event is observed by Geotail at ^^03: 10 UT on 11 
February 1995 at an inbound crossing of the bow shock near 
the subsolar point, i.e. at (12.2, 4.0, 0.6) Re in the GSE co- 
ordinate. Figure 1 shows the overview of main physical pa- 
rameters of the gradual crossing event. The shock transition 
appears as an abrupt change in both the magnetic field data 
(MGF, Kokubun et al., 1994) as well as the plasma data (LEP, 
Mukai et al., 1994). The shock normal direction estimated 
by the Minimum Variance Analysis (Sonnerup and Cahill, 
1967) was (0.94, 0.31, -0.17) consistent with the shock nor- 
mal derived from the semi-empirical bow shock model of 
Peredo et al. (1995). This model is known to give normal 
directions in agreement with those obtained by the timing 
method of multi-spacecraft (Horbury et al., 2002). The up- 
stream parameters. Ma and Obu, were then estimated to be 
^6.8 and ^^68°, respectively. (For various methods of shock 
normal determination, see, e.g. Paschmann and Daly, 1998.) 
Noticeable in the magnetic field data in Figure 1 are the rel- 
atively large fluctuations increasing with time (03:08-03:10 
UT, SB/B ~20%). Our particular interest also goes to the 
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Fig. 1 . From top to bottom are (a) magnitude | B | , (b) azimutlial 4>b , and (c) 
latitudinal 6b component of the magnetic field (MGF) with 3 sec (black 
line) and 1/16 sec (gray line) sampling, and the energy-time (Et) diagrams 
of (d) solar wind ions (LEP/SWI), (e) omni-directional energetic ions 
(LEP/EAI), and (f) omni-directional electrons (LEP/EAE), respectively. 



time profile of higher energy (>0.5 keV) electrons as can be 
seen in the panel (f) as a smooth increase of count rate in the 
time period of 03:09-03: 1 1 UT. It is noted that the spacecraft 
soon exited from the magnetosheath at 03:15 UT. 

Figure 2 shows the detailed spectral properties of the 
waves. They are dominated by the right hand polarized com- 
ponent accompanied by a frequency cut off at around f^H 
which threshold seems not to be pointed out before. We 
speculate that this is because whistler generation concerns 
both ion and electron dynamics. There are tens of detailed 
models for whistler generation and it is not the purpose of 
this letter to discuss the physical meaning of /lh- From the 
Minimum Variance Analysis as well as the Means method 
(Means, 1972), the propagation angle 9kB (the cone angle 
between /c-vector and background magnetic field) were es- 
timated to be 20-40°. We also removed the 180° ambiguity 
of the estimated fc-vector using one component of electric 
field (EFD) data (Matsui et al., 1997). As a result, the waves 
with frequencies lower than ~10 Hz indicated propagation 
toward the sun, away from the shock front, consistent with a 
past report (Orlowski et al., 1994). For higher frequencies, 
we could not obtain reliable results on the propagation direc- 
tion, probably due to the low intensities of the waves. Note 
the cone angle between /c-vector and the solar wind Vsw, 
6kv, were estimated to be 60 - 90° so that the Doppler shift 
was not significant. Supportingly, the spectral slope of the 
high frequency range is approximately ~ 5, consistent with 
past observations of right-hand polarized whistler waves (Or- 
lowski et al., 1995). While traveling upwind, the waves suf- 
fered considerable (exponential) damping as shown in the 
left hand side of Figure 3 which shows the temporal profiles 



of band pass filtered magnetic field data. The characteristic 
time scale of the damping was calculated for each best fit 
model shown by the gray curves. This time scale was found 
to be 47 sec on average. 

All observed features described above are well consistent 
with those of the so-called '1 Hz whistlers' reported else- 
where (e.g. Fairfield, 1974; Sentman et al., 1983; Orlowski 
et al., 1994). From the above arguments, we conclude the ob- 
served waves to be the right hand polarized whistler waves 
propagating away from the shock front. 

It is to be emphasized that, while thermal electrons had 
been studied with respect to whistler wave generation (e.g. 
Tokar et al., 1984; Orlowski et al., 1995), non-thermal elec- 
trons, to the author's knowledge, had never been observed 
in association with the upstream whistlers. Note again that 
the gradual profile of non-thermal electrons, has not been 
analyzed in past in connection with the bow shock crossing 
events as we will discuss below. 

The right hand side of Figure 3 shows the time profiles of 
electron phase space densities (PSDs). The flux increased 
exponentially as the spacecraft approached the shock front. 
The characteristic time scale of the increase was calculated 
for each best fit model shown by the red curves. Above 
2 keV, the typical time scale was 24 sec. During the flux 
increase, a pitch angle distribution also changed as shown 
in Figure 4. Black symbols indicate far upstream region 
of the shock and they show strong asymmetry (i.e. larger 
flux in /I < region, where /i is the cosine of pitch angle 
a), indicating that electrons were streaming away from the 
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Fig. 2. From top to bottom are (a) the magnetic field magnitude (equivalent 
with Figure 1(a)), (b) root mean square of band-pass filtered magnetic 
field, (c) FFT spectra obtained at region A, B, and C indicated by the 
shaded regions in panel (b). The root mean square was obtained from 
frequency above 4 times the ion cyclotron frequency and has been nor- 
malized by the ambient magnetic field magnitude. Arrows in panel (c) 
show ion cyclotron frequency (/ci=0. 19Hz) and lower-hybrid frequency 
(Ilh = •\/ /ci/ce=8. 4Hz), whereas the black curve shows the back- 
ground level obtained from 03:07:05 - 03:07:20 UT (indicated by two 
vertical lines). The red, green, and blue corresponds to the R, L, and C 
component from the fluxgate magnetometer (MGF/FX, <8Hz, fight col- 
ored) and the search coil magnetometer (MGF/SC, <32Hz, heavy col- 
ored), respectively. 
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Fig. 3. (left) Temporal distributions of averaged power spectral densities 
(red, green, and blue for R, L, and C components) and their exponential 
fits (gray lines, 03:08:00 - 03:10:00 UT). Numbers in each panel show 
frequency ranges used for averaging. Frequencies are chosen so that 
they are the correspondence of the resonance conditions of electrons with 
energy shown to the right, (right) Temporal distributions of electrons 
(black lines) and their exponential fits (red lines, 03:08:00 - 03:10:30 
UT). Numbers in each panel show center energies of each energy channel. 



shock front. However, substantial amount of electron counts 
were detected in the /i > region during the time interval 
from 03:08:30 to 03:10:00 UT as indicated by red symbols. 
By this time, the spacecraft was immersed in the precursor 
waves. The distributions isotropized in the immediate down- 
stream (green and blue symbols). 

Figure 5 shows the electron energy spectra. The gray line 
shows the spectrum obtained at 03:00 UT in the pure solar 
wind where there was no contamination from the bow shock. 
The open squares show the spectrum at 03:09:45 UT just 
prior to the crossing. There was a significant amount of en- 
ergetic (>lkeV) electrons compared to the solar wind. The 
spectrum is roughly a power law. The filled squares show the 
spectrum at 03:10:26 UT just within the middle of the shock 
ramp. It now forms a complete power law above 2 keV. Then, 
we applied a chi-square fit above 2 keV with a power-law 
f{E) oc E^^ exp (— SZ-Eioii-off), where E is the electron 
energy as variable, r=4.3(± 0.05), E,.oii-off=3.5(± 0.6), and 
the figures in the parenthesis are the 68% confidence region. 
The resultant /d.o.f. was 28/21, where d.o.f. is the degree 
of freedom and is equal to the number of data points minus 
the number of free parameters. A fit with the kappa distribu- 
tion covering the whole energy range yielded relatively high 
X^/ d.o.f.. Similar spectrum was observed in the immediate 
downstream as well. Note that Eioii-off was introduced to 
better fit the observation but we have not succeeded in deriv- 
ing information of maximum attainable energy of electrons 
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Fig. 4. Pitch angle distributions of electrons of energies, 1.89, 2.25, 2.69, 
3.21, 3.83, and 4.57 keV. The vertical and horizontal axes show PSDs 
and fj, (the cosine of pitch angle a), respectively. Since the magnetic 
field was directed towai'd the sun, fj, > corresponds propagation away 
from the shock front. We have organized the obtained three dimensional 
distribution function by fj, and the horizontal axis has been binned into 16 
bins. The color code indicates different time intervals (in UT) as shown 
in the annotation. 
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Fig. 5. Energy spectra of electrons during the bow shock transition. 



(Oka et al., 2006). 

3. Discussion 

The 'gradual' profile of electron flux, the scattering of 
particles, and the appearance of power-law energy spectrum 
remind us of the diffusive shock acceleration (DSA, e.g. 
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Blandford and Ostriker, 1978). In this process, particles are 
continuously scattered and move back and forth across the 
shock front to gain a net momentum from the velocity dif- 
ference between upstream and downstream. However, there 
are some difficulties in the DSA to fully explain the observed 
features. The observed density ratio r ~ 2.5 {N\ = 21/cm^ 
for the upstream and N2 ~ 52/cm^ for the downstream) and 
the classical DSA formula Tdsa = 3r/2(r - 1) = 2.5 do 
not reproduce the observed T ~ 4.3. The quantitative dif- 
ference can readily be interpreted by, for instance, a free es- 
cape boundary that allows significant number of particles to 
escape from the acceleration region. The discrepancy be- 
tween the e-folding distances derived from wave measure- 
ments (with the resolution of 47 sec) and particle measure- 
ments (with the resolution of 24 sec) indicates that the waves 
were not generated by the non-thermal electrons and that 
self-scattering was weak. Therefore, we do not consider the 
classic DSA mechanism to be fully responsible for genera- 
tion of the observed power-law spectra. 

Nevertheless, we still consider there was substantial scat- 
tering in the upstream as is evident from Figure 4. We have 
shown that the whistler waves were propagating along with 
the majority of the electrons. Since the direction must be 
opposite to that of electron streaming in order to satisfy 
cyclotron resonance condition, the cyclotron resonance by 
the right-hand-polarized whistlers is unlikely as the physi- 
cal mechanism of scattering. The observed whistlers were 
oblique waves (0^3=20-40°), and so the left hand polar- 
ized component and/or electrostatic component of the waves 
might have been playing a role in electron scattering. It is 
also to be mentioned that the amplitudes of the upstream 
waves were relatively high so that a non-linear effect should 
have played a role not only for the scattering but also ac- 
celeration. A recent theory indeed points out a possibility 
of efficient particle acceleration at a turbulent magnetic field 
where aU scales are larger than the particle gyroradius (e.g. 
Giacalone, 2005; Jokipii and Giacalone, 2007). That the 
spectral index was relatively large also favors this theory, al- 
though the detailed comparison is left for future work. 

We focus on the fact that the gradual profile was found at 
the shock with Ma sHghtly below M^^-^ (MaM^^^^ ~0.65 in 
this event). Our interpretation is as follows. In a higher Ma 
shock (Ma >M^^^^), whistler waves do not propagate up- 
stream so that any electron that escapes the shock front can- 
not interact with the waves and that the intense flux can only 
be found at the shock front resulting in 'spike' events. On 
the other hand, in a much lower Ma shock (Ma ^M^;^), 
there would be plenty of waves to scatter particles, but the 
number of non-thermal electrons that are subject to scatter- 
ing is small. Therefore, we expect to find similar events in 
the Mach number range slightly below M^;^. . 

Unfortunately, however, the time resolution needs to be 
sufficiently high to resolve the decaying profile of non- 
thermal electrons. In the 78 events of Oka et al. (2006), 
there were a few events that seemed to be 'gradual', but the 
decaying time was of the order of the time resolution of the 
particle measurement (12 sec), and thus we could not repli- 
cate the analysis presented in this paper. We anticipate that 
more sophisticated observations of multi-spacecraft mission 
such as Cluster, MMS, and SCOPE/CrossScale will reveal 



the nature of 'gradual' events. 
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